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Abstract 

The Klein-Gordon equation in L>-dimensions for a recently proposed Kratzer 
potential plus ring-shaped potential is solved analytically by means of the 
conventional Nikiforov-Uvarov method. The exact energy bound-states and 
the corresponding wave functions of the Klein-Gordon are obtained in the 
presence of the noncentral equal scalar and vector potentials. The results 
obtained in this work are more general and can be reduced to the standard 
forms in three-dimensions given by other works. 
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I. INTRODUCTION 



In various physical applications including those in nuclear physics and high energy physics 
[1,2], one of the interesting problems is to obtain exact solutions of the relativistic equations 
like Klein-Gordon and Dirac equations for mixed vector and scalar potential. The Klein- 
Gordon and Dirac wave equations are frequently used to describe the particle dynamics in 
relativistic quantum mechanics. The Klein-Gordon equation has also been used to under- 
stand the motion of a spin-0 particle in large class of potentials. In recent years, much 
efforts have been paid to solve these relativistic wave equations for various potentials by 
using different methods. These relativistic equations contain two objects: the four-vector 
linear momentum operator and the scalar rest mass. They allow us to introduce two types 
of potential coupling, which are the four-vector potential (V) and the space-time scalar 
potential (S). 

Recently, many authors have worked on solving these equations with physical potentials 
including Morse potential [3] , Hulthen potential [4] , Woods-Saxon potential [5] , Posch-Teller 
potential [6] , reflectionless-type potential [7] , pseudoharmonic oscillator [8] , ring-shaped har- 
monic oscillator [9], Vq tanh^(r/ro) potential [10], five-parameter exponential potential [11], 
Rosen- Morse potential [12], and generahzed symmetrical double- well potential [13], etc. It 
is remarkable that in most works in this area, the scalar and vector potentials are almost 
taken to be equal (i.e., S = V) [2,14]. However, in some few other cases, it is considered 
the case where the scalar potential is greater than the vector potential (in order to guar- 
antee the existence of Klein-Gordon bound states) (i.e., S > V) [15-19]. Nonetheless, such 
physical potentials are very few. The bound-state solutions for the last case is obtained for 
the exponential potential for the s-wave Klein-Gordon equation when the scalar potential is 
greater than the vector potential [15]. 

The study of exact solutions of the nonrelativistic equation for a class of non-central po- 
tentials with a vector potential and a non-central scalar potential is of considerable interest 
in quantum chemistry [20-22]. In recent years, numerous studies [23] have been made in 
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analyzing the bound states of an electron in a Coulomb field with simultaneous presence 
of Aharanov-Bohm (AB) [24] field, and/or a magnetic Dirac monopole [25], and Aharanov- 
Bohm plus oscillator (ABO) systems. In most of these works, the eigenvalues and eigen- 
functions are obtained by means of seperation of variables in spherical or other orthogonal 
curvilinear coordinate systems. The path integral for particles moving in non-central poten- 
tials is evaluated to derive the energy spectrum of this system analytically [26]. In addition, 
the idea of SUSY and shape invariance is also used to obtain exact solutions of such non- 
central but seperable potentials [27,28]. Very recently, the conventional Nikiforov-Uvarov 
(NU) method [29] has been used to give a clear recipe of how to obtain an explicit exact 
bound-states solutions for the energy eigenvalues and their corresponding wave functions in 
terms of orthogonal polynomials for a class of non-central potentials [30] . 

Another type of noncentral potentials is the ring-shaped Kratzer potential, which is a 
combination of a Coulomb potential plus an inverse square potential plus a noncentral angu- 
lar part [31,32]. The Kratzer potential has been used to describe the vibrational- rotational 
motion of isolated diatomic molecules [33] and has a mixed-energy spectrum containing both 
bound and scattering states with bound-states have been widely used in molecular spec- 
troscopy [34]. The ring-shaped Kratzer potential consists of radial and angular-dependent 
potentials and is useful in studying ring-shaped molecules [22]. In taking the relativistic 
effects into account for spin-0 particle in the presence of a class of noncentral potentials, Ya- 
suk et al [35] applied the NU method to solve the Klein-Gordon equation for the noncentral 
Coulombic ring-shaped potential [21] for the case V = S. Further, Berkdemir [36] also used 
the same method to solve the Klein-Gordon equation for the Kratzer-type potential. 

Recently, Chen and Dong [37] proposed a new ring-shaped potential and obtained the 
exact solution of the Schrodinger equation for the Coulomb potential plus this new ring- 
shaped potential which has possible applications to ring-shaped organic molecules like cyclic 
polyenes and benzene. This type of potential used by Chen and Dong [37] appears to be 
very similar to the potential used by Yasuk et al [35]. Moreover, Cheng and Dai [38], 
proposed a new potential consisting from the modified Kratzer's potential [33] plus the 
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new proposed ring-shaped potential in [37]. They have presented the energy eigenvalues 
for this proposed exactly-solvable non-central potential in three dimensional (i.e., D = 3)- 
Schrodinger equation by means of the NU method. The two quantum systems solved by 
Refs [37,38] are closely relevant to each other as they deal with a Coulombic field interaction 
except for a slight change in the angular momentum barrier acts as a repulsive core which 
is for any arbitrary angular momentum £ prevents collapse of the system in any dimensional 
space due to the slight perturbation to the original angular momentum barrier. Very recently, 
we have also applied the NU method to solve the Schrodinger equation in any arbitrary D- 
dimension to this new modified Kratzer-type potential [39,40]. 

The aim of the present paper is to consider the relativistic effects for the spin-0 parti- 
cle in our recent works [39,40]. So we want to present a systematic recipe to solving the 
D-dimensional Klein-Gordon equation for the Kratzer plus the new ring-shaped potential 
proposed in [38] using the simple NU method. This method is based on solving the Klein- 
Gordon equation by reducing it to a generalized hypergeometric equation. 

This work is organized as follows: in section II, we shall present the Klein-Gordon 
equation in spherical coordinates for spin-0 particle in the presence of equal scalar and 
vector noncentral Kratzer plus the new ring-shaped potential and we also separate it into 
radial and angular parts. Section III is devoted to a brief description of the NU method. 
In section IV, we present the exact solutions to the radial and angular parts of the Klein- 
Gordon equation in D-dimensions. Finally, the relevant conclusions are given in section 
V. 

II. THE KLEIN-GORDON EQUATION WITH EQUAL SCALAR AND VECTOR 

POTENTIALS 

In relativistic quantum mechanics, we usually use the Klein-Gordon equation for de- 
scribing a scalar particle, i.e., the spin-0 particle dynamics. The discussion of the relativistic 
behavior of spin-zero particles requires understanding the single particle spectrum and the 
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exact solutions to the Klein Gordon equation which are constructed by using the four- vector 
potential (A = 0, 1, 2, 3) and the scalar potential (S). In order to simplify the solution 
of the Klein-Gordon equation, the four- vector potential can be written as = (^o, 0, 0, 0). 
The first component of the four- vector potential is represented by a vector potential (V), i.e., 
Aq — V. In this case, the motion of a relativistic spin-0 particle in a potential is described by 
the Klein-Gordon equation with the potentials V and S [1]. For the case S >V, there exist 
bound-state (real) solutions for a relativistic spin-zero particle [15-19]. On the other hand, 
for S = V, the Klein-Gordon equation reduces to a Schrodinger-like equation and thereby 
the bound-state solutions are easily obtained by using the well-known methods developed 
in nonrelativistic quantum mechanics [2]. 

The Klein-Gordon equation describing a scalar particle (spin-0 particle) with scalar 
S{r, 9, ip) and vector V(r, 6*, (p) potentials is given by [2,14] 

{p2 _ - V{t, e, ip)/2f +[^1 + S{r, 6, ip)/2f} ^(r, 6, ^) = 0, (1) 

where Er, P and /i are the relativistic energy, momentum operator and rest mass of the 
particle, respectively. The potential terms are scaled in (1) by Alhaidari et al [14] so that 
in the nonrelativistic limit the interaction potential becomes V. 

In this work, we consider the equal scalar and vector potentials case, that is, S{r, 9, (p) — 
V{r, 9, (p) with the recently proposed general non-central potential taken in the form of the 
Kratzer plus ring-shaped potential [38-40]: 



V(.,0,^) = V,M + ^ + i%, (2) 

sm 9 



Vi{r) = + ^, V2{9) = Cctg'9, = 0, (3) 

where A = 2aoro, B = aQr^ and C is positive real constant with gq is the dissociation energy 
and To is the equilibrium internuclear distance [33]. The potentials in Eq. (3) introduced 
by Cheng-Dai [38] reduce to the Kratzer potential in the hmiting case of C = [33]. In fact 
the energy spectrum for this potential can be obtained directly by considering it as special 
case of the general non-central seperable potentials [30] . 



In the relativistic atomic units {h — c — 1), the D-dimensional Klein-Gordon equation 

in (1) becomes [41] 



j.D-1 Qj, \ Qy, 



1 d / . d\ 

^de\ de)^ sin^ e dip^ 



sin^ 9 



with '0(r, 9, ip) being the spherical total wave function separated as follows 



(4) 



Inserting Eqs (3) and (5) into Eq. (4) and using the method of separation of variables, the 
following differential equations are obtained: 



1 d (^u-idR{r) 



r^-^ dr 



dr 



j{j + D-2) 



2/2 ^ ^ 



R{r) = 0, 



(6) 



1 d . ^d' 

I gi Ti H 

sin^d^l d9, 



sin^^ 



+ j{j + D-2) 



H{9) = 0, 



(7) 



(8) 



where = /i — Er, a\= ii + Er, m and j are constants and with and Xj = j{j + D — 2) 
are the separation constants. 

For a nonrelativistic treatment with the same potential, the Schrodinger equation in 
spherical coordinates is 



1 ^ fr^-^D+L 



„D-1 



dr 



dr 



I d ( . 1 

sva.9— + 



sin^a^ 



39 ; sin^ 9 dcp^ 



+ 2/1 



E 



NR 



Viir) 



V2{9) Vsicp) 



sin^ 9 



ip{r,9,^)=0. 



(9) 



where and E^r are the reduced mass and the nonrelativistic energy, respectively. Besides, 
the spherical total wave function appearing in Eq. (9) has the same representation as in Eq. 
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(5) but with the transformation j 
differential equations [39,40]: 



£. Inserting Eq. (5) into Eq. (9) leads to the following 



1 d 



nD-1 



dr 



.idR{r) 
dr 



( A 
- 2/x I^Enr + -- 



R{r) = 0, 



(10) 



1 d / . d_ 

^der^'^de 



+ 2fxC cos^ 9 



+ e{e + D-2) 



H{e) = 0, 



(11) 



d^^if) 
dif^ 



+ m'^^{(p) = 0, 



(12) 



where and Xe = i{i + D — 2) are the separation constants. Equations (6)- (8) have the 
same functional form as Eqs (10)-(12). Therefore, the solution of the Klein-Gordon equation 
can be reduced to the solution of the Schrodinger equation with the appropriate choice of 
parameters: j —>■ i, of —>■ —E^r and — > 2fi. 

The solution of Eq. (8) is well-known periodic and must satisfy the period boundary 
condition -|- 27r) = ^{<p) which is the azimuthal angle solution: 

1 



exp(±im</7), m = 0, 1, 2, 



(13) 



Additionally, Eqs (6) and (7) are radial and polar angle equations and they will be solved 
by using the Nikiforov-Uvarov (NU) method [29] which is given briefly in the following 
section. 



III. NIKIFOROV-UVAROV METHOD 

The NU method is based on reducing the second-order differential equation to a gener- 
alized equation of hypergeometric type [29]. In this sense, the Schrodinger equation, after 
employing an appropriate coordinate transformation s — s{r), transforms to the following 
form: 

"^"^"^ + ^/-^'^ + ^^'^^'^ = °' ^^^^ 
7 



where a{s) and a{s) are polynomials, at most of second-degree, and t(s) is a first-degree 
polynomial. Using a wave function, ipn{s), of the simple ansatz: 

1pn{s) = 0n(s)yn(s), (15) 

reduces (14) into an equation of a hypergeometric type 

ais)y';{s) + T{s)yM + Ay„(s) = 0, (16) 

where 

(17) 

t{s) =r{s) + 2n{s), t'{s) < 0, (18) 

and A is a parameter defined as 

A = A„ = -nr'is) - ~ ^K " (s) , n = 0, 1, 2, .... (19) 

The polynomial t(s) with the parameter s and prime factors show the differentials at first 
degree be negative. It is worthwhile to note that A or A„ are obtained from a particular 
solution of the form y{s) = yn{s) which is a polynomial of degree n. Further, the other 
part yn{s) of the wave function (14) is the hypergeometric-type function whose polynomial 
solutions are given by Rodrigues relation 

where S„ is the normalization constant and the weight function p{s) must satisfy the con- 
dition [29] 

^w{s) = ^M^), w{s) = a{s)p{s). (21) 
as a[s) 

The function tt and the parameter A are defined as 



' '^'W-^M V-g(,) + toM. (22) 



A = A; + 7r'(s). 



(23) 



In principle, since 7r(s) has to be a polynomial of degree at most one, the expression under 
the square root sign in (22) can be arranged to be the square of a polynomial of first degree 
[29]. This is possible only if its discriminant is zero. In this case, an equation for k is 
obtained. After solving this equation, the obtained values of k are substituted in (22). In 
addition, by comparing equations (19) and (23), we obtain the energy eigenvalues. 

IV. EXACT SOLUTIONS OF THE RADIAL AND ANGLE-DEPENDENT 

EQUATIONS 

A. Separating variables of the Klein-Gordon equation 

Wc seek to solving the radial and angular parts of the Klein-Gordon equation given 
by Eqs (6) and (7). Equation (6) involving the radial part can be written simply in the 
following form [39-41]: 



(M-l)(M-3) 
4^2 



a 



A B 



2 2 

+ aia2 



g{r) = 0, 



(24) 



where 



M^D + 2j. 



(25) 



On the other hand, Eq. (7) involving the angular part of Klein-Gordon equation retakes the 
simple form 



^ + ctg{e)^^ 



de^ 



dO 



+ Cal cos^ 9 



-j{j + D-2) 



H{e) = 0. 



(26) 



Thus, Eqs (24) and (26) have to be solved latter through the NU method in the following 
subsections. 
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B. Eigenvalues and eigenfunctions of the angle-dependent equation 



In order to apply NU method [29,30,33,35,36,38-40,42-44], we use a suitable transforma- 
tion variable s = cos9. The polar angle part of the Klein Gordon equation in (26) can be 
written in the following universal associated-Legendre differential equation form [38-40] 



d'^H{s) 2s dH{s) 



+ 



j{j + D -2){1- s') -m' - Cais' H{e) = 0. (27) 



l-s'^ ds ' (l-s2)2 

Equation (27) has already been solved for the three-dimensional Schrodinger equation 
through the NU method in [38]. However, the aim in this subsection is to solve with different 
parameters resulting from the L>-space-dimensions of Klein-Gordon equation. Further, Eq. 
(27) is compared with (14) and the following identifications are obtained 



T(S) 



-2s, a{s) = l-s\ a{s) = -m" + (1 - s')v\ 



(28) 



where 



I/' = + D-2)= j{j + D-2) + Cai = m' + Ca 



2- 



(29) 



Inserting the above expressions into equation (22), one obtains the following function: 



7r(s) = ±^{i^' - k)s^ + k-u' + m'2. 

Following the method, the polynomial 7i{s) is found in the following possible values 

m's for ki — u' — m'^, 
-m's for ki = u' — m'^, 
for k'i — v', 
for k2 — v'. 



(30) 



7r(s) 



(31) 



m 



-m 



Imposing the condition t'{s) < 0, for equation (18), one selects 



ki — u' — m! and 7r(s) = —m's, 



(32) 



which yields 



10 



t(s) = -2(l + m')s. (33) 
Using equations (19) and (23), the following expressions for A are obtained, respectively, 

A = A„ = 2n(l + m')+ri(ri-l), (34) 



A = I/' -m'(l + m'). (35) 
We compare equations (34) and (35), the new angular momentum j values are obtained as 
J = -^-^^Y^ + \sJ{D- 2)2 + (2n + 2m' + 1)2 - ACal - 1, (36) 



or 



/ = + \\I{D- 2)2 + (2n + 2m' + 1)2 - 1. (37) 

Using Eqs (15)-(17) and (20)-(21), the polynomial solution of y„ is expressed in terms of 
Jacobi polynomials [39,40] which are one of the orthogonal polynomials: 

H~{e) = 7V~sin-»Pf ''"'^(cos^), (38) 
where N~ = .^^ /I 2n+2m^+i)(n+2 m0!n! normahzatiou constant determined by 

" 2"* {n+m')\ \ I •' 

+1 2 

/ \B.~{s)\ ds — 1 and using the orthogonality relation of Jacobi polynomials [35,45,46]. 



-1 

Besides 



n _|_ OmM 1 I 

n^- ^ 2 + 2 + + ^^^^ ~ + ^^^^ 



where m' is defined by equation (29). 



C. Eigensolutions of the radial equation 

The solution of the radial part of Klein-Gordon equation, Eq. (24), for the Kratzer's 
potential has already been solved by means of NU- method in [39]. Very recently, using 
the same method, the problem for the non-central potential in (2) has been solved in three 
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dimensions (3D) by Cheng and Dai [36]. However, the aim of this subsection is to solve 
the problem with a different radial separation function g{r) in any arbitrary dimensions. In 
what follows, we present the exact bound-states (real) solution of Eq. (24). Letting 



= afai AY = (M - 1)(M - 3) + ABai (5' = Aai 

and substituting these expressions in equation (24), one gets 

d'^g(r) ( -£^r^ + - Y^ 



(40) 



+ 



9{r) = 0. 



(41) 



To apply the conventional NU-method, equation (41) is compared with (14), resulting in 
the following expressions: 



T(r) — 0, cr(r) — r, (j(r) — —e r + P r — j 
Substituting the above expressions into equation (22) gives 

7r(r) = - ± -^4£2r2 + 4{k - (5'^)r + 472 + 1. 



(42) 



(43) 



Therefore, we can determine the constant k by using the condition that the discriminant of 
the square root is zero, that is 



k^0^± £^472 + 1, 47^ + 1 = (D + 2j - 2f + ABa. 



(44) 



In view of that, we arrive at the following four possible functions of 7r(r) : 



TT r = < 



+ 



+ 



(45) 



er + |V47' + ij for k^ = + e^A^ + 1, 

'er + \yJAY + 1] for A;i = + sVi^TT, 

er - \^AY + 1] for k2 = - e^A^ + 1, 

er - \^JAY + 1] for k2 = (3^ - e^A^ + 1- 

The correct value of 7r(r) is chosen such that the function r(r) given by Eq. (18) will have 
negative derivative [29]. So we can select the physical values to be 



k = /3- eJAY + and 7r(r) 



(46) 
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which yield 



T(r) = -2er + (1 + ^472 + 1), r'(r) = -2e < 0. (47) 
Using Eqs (19) and (23), the following expressions for A are obtained, respectively. 



A = A„ = 2ne, n = 0, 1, 2, (48) 



A = 5'-£(l + V47' + l). 



(49) 



So we can obtain the Klein Gordon energy eigenvalues from the following relation: 



1 + 2n + V(^ + 2j - 2y + 4(/x + Er)B 



H - Er = AJ 11 + Er, 



and hence for the Kratzer plus the new ring-shaped potential, it becomes 



l + 2n+J{D + 2j - 2Y + iaoriii^ + Er) 



jji- Er = 2aoroJiJ, + Er, 



(50) 



(51) 



with j defined in (36). Although Eq. (51) is exactly solvable for Er but it looks to be little 
complicated. Further, it is interesting to investigate the solution for the Coulomb potential. 
Therefore, applying the following transformations: A — Ze^, B — 0, and j — i, the central 
part of the potential in (3) turns into the Coulomb potential with Klein Gordon solution for 
the energy spectra given by 

^ 2q^e^ ^ 



Er^I^II 



n,£ = 0,l,2,.... 



(52) 



^ g2g2 + (2n + 2£ + L>-l)2^ 

where q — Ze is the charge of the nucleus. Further, Eq. (52) can be expanded as a series in 
the nucleus charge as 



Er^ II- 



+ 



-3W - 0{qe)\ 



(53) 



(2n + 2£ + L'- 1)2 {2n + 2e + D 

The physical meaning of each term in the last equation was given in detail by Ref. [36]. 

Besides, the difference from the conventional nonrelativistic form is because of the choice of 
the vector V{r, 9, (f) and scalar S{r, 9, (f) parts of the potential in Eq. (1). 
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Overmore, if the value of j obtained by Eq.(36) is inserted into the eigenvalues of the 
radial part of the Klein Gordon equation with the noncentral potential given by Eq. (51), 
we finally find the energy eigenvalues for a bound electron in the presence of a noncentral 
potential by Eq. (2) as 

1 + 2n + ^{2f + D-2f + 4(aor2 - C)(// + Er)] ^h-Er = ^a^r^^ n + Er, (54) 



where m' — \J m?' + C(// + Er) and n is given by Eq. (39). On the other hand, the solution 
of the Schrodinger equation, Eq. (9), for this potential has already been obtained by using 
the same method in Ref. [39] and it is in the Coulombic-like form: 

Enr = n = 0, 1, 2, ... (55) 



2n + 1 + ^(2f + L> - 2)2 + 8At(aorg - C) 



2i' + D-2 = ^{D - 2)2 + (2n + 2m' + ly - 1, (56) 

where m' = \Jm^ + 2\xC. Also, applying the following appropriate transformation: ij,+Er — >■ 
2/1, II — Er — > — Er^R, j ^ £ to Eq. (54) provides exactly the nonrelativistic limit given by 
Eq. (55). 

In what follows, let us now turn attention to find the radial wavefunctions for this 
potential. Substituting the values of a{r),n{r) and T(r) in Eqs (42), (45) and (47) into Eqs. 
(17) and (21), we find 

0(r) = r(^+^)/2g-er^ (57) 

p{r) = r^e-2^^ (58) 

where ( = ^47^ + 1. Then from equation (20), we obtain 

y^,{r) = i?„,r-^e^-— (r^+^e-^-) , (59) 

and the wave function g{r) can be written in the form of the generalized Laguerre polyno- 
mials as 
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(1+0/2 



(60) 



where for Kratzer's potential we have 



C = \/{D + 2j - 2f + 4aor2(/x + Er), p = 2er. 



(61) 



Finally, the radial wave functions of the Klein-Gordon equation are obtained 



R{p) = C, 



nj 



2eJ 



(C+2-£))/2 



e-'^'Liip), 



(62) 



where C^j is the normalization constant to be determined below. Using the normalization 

oo 

condition, / R^{r)r^~^dr = 1, and the orthogonality relation of the generalized Laguerre 



polynomials, Z^^+^e"^ [Ll{z)f dz = we have 



C,,^ = [2d p? - El 



V (2n + C + l) (n + C)!' 
Finally, we may express the normalized total wave functions as 



1+1 



(2n + 2m' + l){n + 2m'y.n\n\ 
2'^'{n + m')\ \ 27r(2n + C + l)(n + C)! 



(63) 



xr 



(C+2-D) 



— exp{-Jp^ - Elr)Li{2Jp^ - E^r) sin'">)P^™'''"')(cos ^) exp{±im^). (64) 



where ( is defined in Eq. (61) and m' is given after the Eq. (54). 



V. CONCLUSIONS 

The relativistic spin-0 particle £)-dimensional Klein-Gordon equation has been solved 
easily for its exact bound-states with equal scalar and vector ring-shaped Kratzer potential 
through the conventional NU method. The analytical expressions for the total energy levels 
and eigenfunctions of this system can be reduced to their conventional three-dimensional 
space form upon setting D = 3. Further, the noncentral potentials treated in [30] can 
be introduced as perturbation to the Kratzer's potential by adjusting the strength of the 

15 



coupling constant C in terms of oq, which is the couphng constant of the Kratzer's potential. 
Additionally, the radial and polar angle wave functions of Klein-Gordon equation are found 
in terms of Laguerre and Jacobi polynomials, respectively. The method presented in this 
paper is general and worth extending to the solution of other interaction problems. This 
method is very simple and useful in solving other complicated systems analytically without 
given a restiction conditions on the solution of some quantum systems as the case in the 
other models. We have seen that for the nonrelativistic model, the exact energy spectra can 
be obtained either by solving the Schrodinger equation in (9) (cf. Ref. [39] or Eq. (55)) or by 
applying appropriate transformation to the relativistic solution given by Eq. (54). Finally, 
we point out that these exact results obtained for this new proposed form of the potential 
(2) may have some interesting applications in the study of different quantum mechanical 
systems, atomic and molecular physics. 
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